Evolution of the stripe phase as a function of doping 
from a theoretical analysis of angle-resolved photoemission data 
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By comparing single-particle spectral functions of t-J and Hubbard models with recent angle- 
resolved photoemission (ARPES) results for LSCO and Nd-LSCO, we can decide where holes go 
as a function of doping, and more specifically, which type of stripe (bond-, site-centered) is present 
in these materials at a given doping. For dopings greater than about 12% our calculation shows 
furthermore that the holes prefer to proliferate out of the metallic stripes into the neighboring anti- 
ferromagnetic domains. The spectra were calculated by a cluster perturbation technique, for which 
we present an alternative formulation. Implications for the theory for high-Tc superconductivity 
are discussed. 

PACS numbers: 74.72.-h, 79.60.-i, 71.27.-fa 
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I. INTRODUCTION 



At present, stripes are at the heart of the debate con- 
cerning the mechanism of superconductivity in high- 
temperature superconductors (HTSC). There is clear ex- 
perimental evidence for static stripes in Nd-doped LSCO 
(Nd-LSCO, for example Lai.48Ndo.4Sro.i2Cu04) from 
elastic neutron scattering experiments |0] . The existence 
of a dynamical stripe phase in the "real" high-Tc com- 
pound LSCO or even in YBCO has been conjectured 
from similar diffraction patterns in the inelastic neu- 
tron scattering results (2||^. However, it has not been 
decided so far, whether the observed pattern is due to 
one-dimensional spin-inhomogeneities (i.e. stripes) B] 
or to two-dimensional incommensurable spin-waves M. 
In the case of YBCO near optimal doping, it has been 
suggested that the incommensurable fourfold neutron 
scattering peak is due to the dispersion of the famous 
41meV commensurable (found below Tc at momentum 
k = (tt, tt)) neutron scattering peak to lower excitation 
energies around k = (tt, tt) Q. From the theoretical 
point of view, several numerical analyses, ranging from 
Hartree-Fock ||,§, DMRG to dynamical mean-field 
theory ^ indicate that stripes can be produced by purely 
strong-correlation effects. On the other hand, structural 
transformations [p",pHlO|] as well as long-range Coulomb 
interactions [ pT| , p^ may also play an important role in 
the formation of stripes. 

In this article, we provide numerical arguments show- 
ing that an essential link in the chain of evidence for 
stripes is provided by angle-resolved photoemission spec- 
troscopy (ARPES): ARPES spectra show hardly any dif- 
ference between LSCO (c??/namica/ stripes candidate) and 
Nd-LSCO {static stripe system) [|l§-|l5|. This fact was 
first pointed out by a semiphenomenological argument by 
Salkola et al. p^ . In a previous paper, we have shown 
that salient spectral features of Nd-LSCO and LSCO can 
be explained by a model with static stripes Here, 
we will show that the spectra of LSCO and Nd-LSCO 



can be almost quantitatively described by different types 
of stripe states (i.e. site-centered, bond-centered) for a 
wide variety of dopings. Since the experimental ARPES 
spectra for LSCO and Nd-LSCO are so similar and their 
spectra can be quantitatively described by stripe mod- 
els, we argue that there must be stripes present in LSCO 
as well, at least in the underdoped to optimally doped 
regime. 

If stripes are present as low-energy excitations in the 
high-Tc compounds, they must affect the microscopic de- 
scription of the superconducting state, independently on 
whether they are an obstacle against superconductivity 
or even its driving force. In this context, an important 
question is, whether stripes are bond-centered or site- 
centered since, theoretically, bond-centered stripes have 
been shown to enhance superconducting pairing correla- 
tions | |l8| , p^ . By analyzing the different scenarios with 
our technique and comparing the results with ARPES 
spectra for different dopings, we can decide which type of 
stripe (bond- or site-centered) better describes the spec- 
trum as a function of doping. 

Very recently, Zhou et al. succeeded in measuring 
ARPES data on Nd-LSCO and LSCO systems with dif- 
ferent doping levels. In the doping region beyond 12% 
the ARPES results reveal a dual nature of the electronic 
structure: The straight segments forming the "Fermi sur- 
face" in energy-integrated photoemission and the (tt, 0) 
low-energy excitations, which have been attributed to 
site-centered stripes for the 12% doping case ||l^, are 
still present but overlaid with more two-dimensional 
(2D) features reminiscent of a simple tight-binding band- 
structure for homogeneous 2D systems. Zhou et al. ad- 
dress the experimental question, whether the two fea- 
tures originate from the mixing of two different phases or 
whether they are intrinsic properties of the same stripe 
phase. In the first case (phase separation) there should a 
non-stripe phase with a very high carrier concentration 
(as estimated from the area inside the Fermi surface) 
coexisting with a phase of site-centered stripes. The au- 
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thors remark that there is hardly evidence for such a 
phase in Nd-LSCO and LSCO at the doping levels un- 
der consideration Another scenario consistent 
with the experimental observation would be that the sys- 
tem forms more and more bond-centered stripes upon 
increasing the doping beyond 12%. The possibility for 
the latter scenario was already indicated by our earlier 
calculations where the spectrum of bond-centered 
stripes resembles the diamond-shaped two-dimensional 
feature observed by Zhou et al.. In the present article we 
will show that the evolution of the experimental ARPES 
with increasing doping can be described by assuming that 
more and more bond-centered stripes are formed at the 
expense of site-centered ones (as conjectured previously 
by Zhou et al.). 

Our approach adopts the Cluster-Perturbation-Thcory 
(CPT) developed by Senechal et al. which consists 
in splitting the infinite lattice into clusters which are 
treated by exact diagonalization. The inter-cluster hop- 
ping terms are then treated perturbatively, so that one 
eventually approximates the infinite lattice. This method 
takes into account exactly the local correlations, which 
probably are the most important ones in these systems, 
and at the same time makes all k points of the Brillouin 
Zone available. In addition, this is an ideal method to 
deal with a "larger unit cell" such as the one present in 
the stripe phase. One should, however, mention that, 
at this level, the method is not appropriate to explore 
stripe stability for a given model. As a matter of fact, 
this is not the aim of the present work. Rather, via this 
CPT approach we enforce a stripe pattern by connect- 
ing clusters with different hole dopings [see Fig. (^], and 
study its spectral function to compare with ARPES ex- 
periments. A similar approach has been taken in Ref. 

where the hole spectral function for site-centered 
stripe patterns was calculated within the string picture. 
Our paper is organized as follows: In section ||, we ex- 
plain in detail how the CPT is applied to the stripe phase. 
In section [II our numerical results are presented and 
compared with ARPES spectra. The content of this ar- 
ticle is summarized in section EV. Finally, the appendix 
presents an alternative derivation of the CPT equation, 
by means of a mapping onto a hard-core fermion model. 
In this appendix, we show that the CPT approximation 
amounts to neglecting two-particle excitations within this 
model. 



II. TECHNIQUE 

The computational technique for our calculation of 
the single-particle spectral weight A(k, lo) and for the 
Green's function is a special application of the CPT for 
inhomogeneous systems. This method is based on a 
strong-coupling perturbation expansion of the (Hubbard 



model's) one-body hopping operators linking the individ- 
ual unit-cells At lowest order in this expansion, the 
Green's function of the infinite lattice can be expressed 
(in matrix form) as: 

Qcluster (^^^-^ 



G°°(P,z) = 



1 



(1) 



where the matrices G°°, G'^'"**'^'', and e(P) will be defined 
around Eq. (|^). The authors of ref. have pointed 
out that the above formula becomes exact in the limit 
of vanishing interaction {U /t = 0) and, obviously, in the 
atomic limit {t = 0), and can thus be considered as an 
interpolation scheme between t oo and f — > 0. No- 
tice, however, that this formula does not become exact 
in the C/ — > oo limit. In the interesting regime where 
the interaction and the hopping are of the same order of 
magnitude Senechal et al. have shown numerically that 
Eq. (|l]) gives as an accurate interpolation between the two 
limiting cases. Additionally, we have shown in our pre- 
vious article that the cluster perturbation technique 
is ideally suited to study inhomogeneous systems such as 
the stripe state in the high-Tc compounds. 
In order to deal with stripes, the infinite lattice is di- 
vided into unit cells of equal size (Fig.(|l|a)). The unit 
cell is further divided into independent blocks to incor- 
porate the stripe topology: In the example of Fig.(|^), 
we are interested in a site-centered ("3+1") configura- 
tion with quarter-filled metallic chains alternating with 
half-filled 3-leg ladders. The 3-leg ladders (here 3x6 
blocks) and quarter-filled chains (here 1 x 12 with 6 holes, 
see Fig.(|l|D)) are solved by exact diagonalization, yield- 
ing the single-particle spectral function of the block. The 
individual Green's functions of the blocks forming a unit- 
cell cluster are combined to form the Green's function of 
the unit cell G'^'"^*'^'' at "order zero", i. e., in which the 
intra-cluster hopping terms are set to zero. In a second 
step, the intra-cluster hopping connecting the individual 
blocks both within the same unit cell and in different 
unit cells (dashed lines in Fig.(|^b)) are incorporated via 
the cluster perturbation technique. This yields the de- 
sired Green's function of the infinite lattice. Where it 
was technically feasible, we doubled the unit cell (as in 
Fig.(|l|)) and diagonalized two 3-leg ladders with a stag- 
gered magnetic field pointing in opposite directions, re- 
sulting in a TT-phase shifted (between the AF domains) 
Neel order in the final configuration. This site-centered 
"3 -I- 1" configuration with 7r-phase shifted Neel order of 
stripes was first suggested by Tranquada et al. ||l| . Bond- 
centered stripes, on the other hand, arc modeled by 2-leg 
ladders with alternating filling. In the following, we will 
refer to this bond-centered configuration as "2 -f 2" . 
Holes can propagate out of the metallic stripes into the 
AF insulating domains via the inter- and intra-unit-cell 
hoppings (dashed lines in Fig.(|l|)). As explained above, 
our method consists in "forcing" the stripe structure "by 
hand" in order to study the effects of this structure on 
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the photoemission spectrum. The different hole concen- 
tration in the "metallic" and in the "antiferromagnetic" 
region is achieved by adjusting independently the chem- 
ical potential of the individual blocks so that the desired 
hole density (see Fig. ^) is obtained. This corresponds 
to introducing an on-site energy shift between the two 
blocks A « 1.5t. Physically, this shift corresponds to the 
energy of the stripe formation, possibly produced by a 
combined effect of strong correlation and of lattice distor- 
tion occurring in the low-temperature tetragonal (LTT) 
phase Of course, without such an energy shift it 
would be impossible to obtain such large density oscil- 
lations, such as the ones shown in Fig. [J. Indeed, in a 
homogeneous system the amplitude of charge oscillations 
remains of the order or less than 0.1, as shown, e. g., by 
DMRG calculations 

Alternatively equation(nf) can be obtained by expressing 
the fermionic creation/annihilation operator c|^^ in terms 
of fermionic creation and annihilation operators da rep- 
resenting the photoemission and inverse photoemission 
target statesja) of the diagonalized cluster. This deriva- 
tion of Eq. (|l|) is presented in the appendix. 
In Eq. (|l]) , P is a superlattice wave vector and G°° is the 
Green's function of the "cx)-size" 2D system, however, 
still in a hybrid representation: real space within a cluster 
and Fourier-space between the clusters. This is related 
to the fact that G°°(P, z) is now an M x M matrix in the 
space of site indices (in the inhomogeneous stripe configu- 
ration of Fig. (0b) M = 2x(3x6)-hlxl2= 96). Likewise, 
e(P) and Gciuster ^^.^ M x M matrices in real space with 
e(P) standing for the perturbation. For the situation 
in Fig.(^), the only nonzero elements of the hermitian 
matrix e(P) are represented by the dashed bonds in the 
figure: 

{— t : dashed bonds inside cell 
— te^''^="/« : bonds connecting cells 
: elsewhere 

(2) 

In order to facilitate diagonalization of the individual 
clusters, we used periodic boundary conditions along the 
stripe direction. In principle, this introduces hopping 
terms, which are not present in the infinite lattice. How- 
ever, this is not a problem, since these terms can be 
consistently removed perturbatively by subtracting cor- 
responding terms to the matrix elements of e(P). 
A complete Fourier representation of G°° in terms of the 
original reciprocal lattice then yields the cluster pertur- 
bation theory (CPT) approximation |po| . 
To allow for larger block sizes, we actually diagonalized 
the t-J model on the blocks to obtain the (block-) Green's 
functions as an approximation for the Hubbard model's 
(block-) Green's function. The t-J Hamiltonian is de- 
fined as 



H^-t (4A.+1^-C-) + JE(S»S,"^). (3) 
{i,i>>o" (»>i> 

The sums run over all nearest neighbor pairs No 
double occupancy is allowed. We have chosen the com- 
monly accepted values J/t — 0.4 and t w 0.5eV. 
The quality of this approximation is tested by comparing 
single-particle spectral functions of 2-1-2 bond-centered 
stripe configurations at 12% doping based on diagonal- 
izations of both Hubbard- and t-J models. In Fig.(||), 
we show A{k, oj) for the standard walk through the Bril- 
louin zone. One observes that the result for the Hub- 
bard model {U = 8t, Fig. (j^a)) is very similar to the 
t-J model result (J — QAt, Fig.^D)): The dispersion is 
two-dimensional, metallic-like and comparable to a tight- 
binding dispersion. The only difference between the fig- 
ures is quantitative, namely, the different Fermi velocities 
and bandwidths. This difference is predominantly due to 
the omission of the conditional hopping terms ^ J |22| in 
the t-J Hamiltonian. Another reason may be due to the 
fact that the basic parameters of the models, the interac- 
tion strengths J and U have not been fine-tuned to match 
each other. The integrated spectral weight obtained from 
the two models is depicted in Fig.(||). In both cases, the 
occupation in momentum space n{k) is distributed al- 
most isotropically around the F point. The low-energy 
excitations for both models indicate a two-dimensional 
LDA-like Fermi surface, however with additional weight 
at the (±7r,0), (0, ±7r) points. The additional (7r,0) fea- 
tures are sharper in the Hubbard-modcl case. 
Turning on the perturbation allows the holes to travel be- 
tween the blocks and unit-cells. In Fig.(||), we show the 
electron concentration (averaged in the direction along 
the stripes) in the direction perpendicular to the stripes 
before (thick lines) and after (bars) applying CPT for 
the stripe configurations that are discussed in this arti- 
cle. As can be seen, the holes do not travel far from the 
domains that were originally defined, and the electron 
occupation hardly changes from the unperturbed setup. 
Therefore, the desired stripe configuration is conserved 
in our approach. 



III. NUMERICAL RESULTS 

We proceed to the discussion of the spectra for the under- 
doped region: At 10 % doping, the stripes have a charge 
periodicity of 5 lattice constants according to the Tran- 
quada picture. This configuration has been modeled by a 
5x12 unit cell consisting of two half-filled 4x6 systems 
(stacked on top of each other) and a quarter-filled 12- 
site chain as displayed in Fig.(^). Fig.(Ha) shows that 
the spectral weight is confined in one-dimensional seg- 
ments of the Brillouin zone indicating a one-dimensional 
Fermi surface. The low-energy excitations in Fig. (Qa) are 
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mainly located at the (±7r, 0) and (0, ±7r) points in mo- 
mentum space. In Fig.(^) the single-particle spectral 
function for this doping is plotted directly and one ob- 
serves the characteristic stripe features that have been 
discussed in detail in a dispersionless band near 

(tt, 0) crossing the Fermi surface at (tt, 7r/4) resulting 
from the one-dimensional chain oriented in y-direction, 
the double-peak structure at (tt, 0) from the hybridiza- 
tion of the metallic band with the top of the antiferro- 
magnetic band and an excitation at (7r/2,7r/2) at higher 
binding energies than at (tt, 0). The combined results for 
this doping agree very well with the recent ARPES re- 
sults by Zhou et al. jlj] giving further support for the 
static stripe picture by Tranquada: at least below 12% 
doping, charge carriers are only present in quarter-filled 
chains that are alternating with half-filled antiferromag- 
netic domains and an increasing of doping is realized by 
lowering the distance between the chains and therefore 
reducing the effective size of the antiferromagnetic do- 
mains. 

The incommensuration of the quasi-elastic Neutron scat- 
tering does not increase any more beyond a doping level 
of 12%. The simple picture of one-dimensional chains 
moving closer together at the expense of antiferromag- 
netic undoped domains thus cannot be valid in this 
regime. For a doping of 1/8, we have previously shown 
I prt that the ARPES data can be explained by assuming 
that the static stripe system (Nd-LSCO) is in a state of 
site-centered stripes, where 3-leg antiferromagnetic lad- 
ders alternate with quarter- filled chains. Here we will 
address the question how to describe the system for dop- 
ings higher than 1/8 or 12%. 

The ARPES results by Zhou et al. suggest that LSCO 
and Nd-LSCO samples at 15 % doping are in a state that 
is still mainly in a 3-f 1 stripe phase since the ARPES re- 
sults show all the features that have been previously ob- 
served for the X = 0.12 samples: one-dimensional Fermi 
surface and low energy excitations located at (±7r, 0) 
and (0,±7r). In addition, however, low-energy excita- 
tions appear around the edge of a diamond located at 
the center of the Brillouin zone. These excitations are 
connecting the (tt, 0) features. Zhou et al. conjecture 
that bond-centered stripes are formed at the expense of 
site-centered ones, since it was shown in our previous 
calculation that bond-centered stripes do indeed exhibit 
such a diamond-shaped low-energy excitation pattern. 
This description is particularly interesting since bond- 
centered stripes in contrast to site-centered ones have 
been shown to enhance superconducting pairing corre- 
lations p8[ . Therefore, this picture may provide a link 
to the doping dependence of the superconducting transi- 
tion temperature of LSCO. However, to be able to relate 
bond-centered stripes to the diamond-shaped low-energy 
excitations that appear in LSCO and Nd-LSCO at 15 % 
doping, one has to study bond-centered stripes at higher 
dopings than x = 15%, since, in a phase-separated state. 



they have to carry all the additional holes that make up 
the overall doping of x = 15% of the experimental sam- 
ple. Here, we study three possibilities of "overdoped" 
stripe configurations as displayed in Fig.(^ (labeled b,c,d 
consistent with the figure label of Ref. |lj]): 

b) 2+2 bond-centered configuration (AF doped), 
Fig.(^): Here, one ladder is at the filling (n) = 0.75 as 
in the case of 12% doping and the other ladder (previ- 
ously undoped in the 12% doping case) with a filling of 
{n) = 0.875 yielding an overall doping of 19 %. In this 
scenario the doped region extends into the antiferromag- 
netic domain between the charged stripes. Technically, 
this configuration has been realized by coupling a 2 x 8 
ladder (2 holes) with a 2 x 8 ladder (4 holes). 

c) 2+2 bond-centered configuration, Fig.(^): 
Here, one ladder is at the filling {n) — 0.625 and the other 
ladder stays half-filled as in the 12% doping case, again 
yielding an overall doping of 19 %. In this scenario the 
excess holes further populate the charged stripe. Techni- 
cally, this configuration has been realized by coupling a 
2x8 ladder (no holes) with a 2 x 8 ladder (6 holes). 

d) 3+1 site-centered configuration, Fig.(^): Here, 
the chain is quarter-filled ((n) — 0.5) as in the case of 12% 
doping and the three-leg ladder (previously undoped in 
the 12% doping case) has a filling of (n) ~ 0.89 yielding 
an overall doping of 21 %. In this scenario, the doped re- 
gion extends into the antiferromagnetic domain between 
the charged stripes (as in case b). Technically, this con- 
figuration has been realized by coupling two 3x6 ladders 
(2 holes each) with a 1 x 12 chain (6 holes). 

In Fig.(^ b,c,d) the electron occupation in momentum 
space is displayed. Neither of the three configuration 
shows the typical stripe signatures (one-dimensional dis- 
tribution of spectral weight) but the weight is more or less 
isotropically distributed around the F-point. The spec- 
tral weight in Fig.(^ c) is almost circularly distributed 
and resembles n{k) of the free electron gas. The weight 
distribution in Figs. (^b,d), on the other hand, is more of 
a diamond shape. In the "phase-separated" picture, we 
expect the spectral weight stemming from the above con- 
figurations (b,c,d) to be superimposed onto the spectral 
weight coming from the dominating 3+1 site-centered 
structure which is all concentrated in the one-dimensional 
segments in momentum space. This might be the reason 
that in the actual experiment Zhou and coworkers can 
only resolve the latter. 

For the low-energy excitations, the situation is differ- 
ent since a much smaller energy window is integrated 
and therefore the experiment is more sensitive to smaller 
amounts of spectral weight: Here, the results of the 
calculation are shown in Figs. (|^ b,c,d). The low- 
energy excitations of the 2+2 configuration c are concen- 
trated around the (±7r/2, ±7r/2) points. Overlaying these 
excitations with the 3+1 site-centered features would 
not yield the experimentally observed diamond structure 
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connecting the (tt, 0) excitations and therefore this con- 
figuration can be discarded. In contrast, the other 2+2 
configuration b, where the extra holes are populating the 
antiferromagnetic domains does indeed show the exper- 
imentally observed diamond shape. In this setup, the 
(tt, 0) features are present as well. They further enhance 
the low-energy excitations in this region of momentum 
space which are due to the domains that are still in the 
3-1-1 site-centered configuration (present at 12% doping). 
The doped 3-1-1 site-centered configuration d, where the 
holes extend into the antiferromagnetic region also has 
its low-energy excitations distributed around a diamond 
centered at the F point (Fig. (0d). However, its features 
are not as sharp as in Fig.(0b). This configuration can- 
not be so easily discarded and might be present in the 
actual material. 

For completeness, we display the single-particle spectral 
weight of the three "overdoped" stripe configurations in 
Figs. (^ b,c,d): Configurations b and d show a two- 
dimensional tight-binding-like dispersion, the difference 
being the large amount of spectral weight that is concen- 
trated near the Fermi level at (tt, 0) for the bond-centered 
configuration b, which was also visible in the low-energy 
excitation plot (Fig.(^ b)). The other bond-centered con- 
figuration c, where all doped holes are concentrated in 
one ladder, only has a Fermi-level crossing at (7r/2, 7r/2) 
("hole-pockets") consistent with Fig.(0 c). 

IV. CONCLUSION 

In this work, the cluster perturbation technique and its 
application to the stripe phase of high-Tc materials and 
related compounds has been studied in detail. The CPT 
has been applied to Hubbard and t-J systems with differ- 
ent dopings and stripe configurations. The comparison 
of our results with recent ARPES data suggest that, in 
the case of LSCO and Nd-LSCO, stripes are present over 
a wide doping range. For dopings below 12% the system 
consists of site-centered stripes, whereas for higher dop- 
ings more and more bond-centered stripes are present at 
the expense of site-centered ones. In the case of bond- 
centered stripes at dopings higher than 12%, we have 
provided evidence that the excess holes prefer to prolif- 
erate out of the stripes into the AF domain. 
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APPENDIX 

The CPT expression for Green's function of the infinite 
lattice, Eq. (P, can be obtained in an alternative way 
by mapping the fermionic operators onto "quasiparticle" 
operators which create the ionization and affinity states 
from the ground state |23y24| . We discuss here for sim- 
plicity the homogeneous case, whereby the lattice is di- 
vided into equal clusters, although the extension to our 
case of different clusters is straightforward. Consider now 
the ground state | G) of the cluster Hamiltonian (say with 
N particles) and the excited states \a) with -f 1 and 
A^ — 1 particles with corresponding excitation energies 
(including chemical potential) We can think the \a) 
as being created from the ground state via a fermionic 
creation operator 

N=4|G), (4) 
which must satisfy the hard-core constraint 

Y,did^<\. (5) 

a 

If one neglects two-particle excitations, the original 
fermion operators cj^ (n is the combined site and spin 
index within the cluster) can be expressed in terms of 
the djj and as 

with 

Tn,ti = {G\Cn\^l) 

Sn,u = {y\Cn\G) . (7) 

Here and in the following, we are using the label /i for 
N + 1-particle (inverse photoemission) states, and v for 
N ~ 1-particle (photoemission) states. Other labels will 
not distinguish between them. In terms of the Cn, the 
inter-cluster hopping part of the Hamiltonian has the 
general form 

T = Y.Y.^n,n{a,h)ci{a)c,^{h) (8) 

n,m a, 6 

where a, h label the individual clusters. Neglecting the 
two-particle excitations and the constraint, the Hamilto- 
nian of the infinite lattice becomes 

H^H, + f, (9) 

where is the intra-cluster Hamiltonian 

7?, = ^^£„4(a)4(a), (10) 
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and we have introduced a cluster index a for the da ■ The 
Hamiltonian Eq. is now quadratic in the da operators 
and can be readily solved exactly by Fourier transforma- 
tion of the inter-cluster part T in the Cluster variables 
a, b. 

We now show that the resulting Green's function for the 
c„ operators is given by Eq.(|l]). For convenience, we first 
carry out a particle- hole transformation on the operators 

4, 



(11) 



such that the new operators all create particles. Eq. 
(O) simplifies to 



~ y ] Qn,a Pq I 



(12) 



where the matrix Qn.a is given in terms of the T and S 
matrices in Eq. (@) as Qn,tj. = T'n,^ and Qn,u = Sn,v ■ The 
particle-hole transformation affects He which becomes 



He = '^£a'na'^Pa{a)Pa{0') + COnst. 



(13) 



with 



^M = +l 
jyi' = -1 : 

and T is transformed to 



(14) 



T' = ^^^%i.m{a,b) Ql^aQ,n,fjPi{o)Pf3{b). (15) 
a,b n,m a,l3 

The total Hamiltonian Eq. (^ can thus be written in the 
form 



H = ^^pl{a)ha,i3{a, b)pp{h) 

a^b a.,f3 



(16) 



with 



ha,p{a, b) = Sa,pSar]a + {Q''T{a, b)Q)^ ^ . (17) 
The Green's function for the p is readily evaluated from 

Eq. m 



{{pa{a),pl{b))) = [z- h)^lji^ , 



(18) 



where we have considered the terms within braces as ma- 
trices in the indices aa, f3b, whereby the complex fre- 
quency z is proportional to the identity matrix. The 
Green's function for the "true" particles c is readily ob- 
tained by inserting the transformation Eq. ( [T^ ) 



GTn,bm ^ (Ma), elm = (Q (^ - h)-' gt) 

V / aa. 



/3b 



(19) 



We now want to show that Eq. ( pj| ) is equivalent to 
Eq. §). 

Consider first the exact Green's function of a single Clus- 
ter 



G'l",fr^'^a,6((c„(a),4(a))),. 



(20) 



By inserting Eq. (|^) in its Lehmann representation, one 
obtains 



G 



cluster 



(21) 



where the matrices e and 77 are diagonal matrices in the 
indices a,/? with values and rja, respectively. 
Using the anticommutation rules, it is now straightfor- 
ward to show that the product Q is equal to the 
identity matrix /. Let us assume for a moment that Q 
is a square matrix, so that Q is a unitary matrix, and 
Q^Q — I as well. In this way, Eq. (19) can be readily in- 
verted yielding (we consider the matrix Q as the identity 
matrix in the a, b indices, i. e. Qna,pb = 5a,bQn,i3) 

^Q{z-h)Q^ 

= Q{z ^erj- Q''TQ)Q'' = G"^''''"'' ~^ - T (22) 

which is equivalent to Eq. (|^), if one considers that e(P) 
(Eq. ^) is nothing but the Fourier transform of T in the 
cluster coordinates a and b. 

Unfortunately, Q is in general not a square matrix, as 
there are more single-particle excited states as particles. 
Therefore, one cannot easily invert Eq. (pj|). Neverthe- 
less, this problem can be readily overcome. We sketch 
the main procedure below. The matrix Qn,a, say with 
dimensions Rx S, consists of R row orthonormal vectors 
f . By appending the remaining S—R orthonormal vec- 
tors to Q, one obtains a square matrix Q which is now 
unitary. For the "extended" Green's functions G°° and 
Qciuster^ obtained by replacing Q ^ Q in Eqs. ( p^ and 
(pT|), respectively, the relation Eq. ( p2| ) obviously holds. 
It is now a matter of matrix algebra to show that G°° and 
Qciuster g^j.g given by the "upper left" R x R blocks (in 
the n, m indices) of the respective "extended matrices, i. 



,bm 



Go 
n 



for n,m < R . 



(23) 



and the same for G"'"'*'"'' . The last line of Eq. (|2|), thus, 
holds for the case of non-square matrices Q as well. 
In summary, we have shown that the CPT is equivalent to 
a mapping onto a model of hard-core fermions describing 
single-particle excitations from the ground state of the 
cluster. This fact suggest an improvement of the method 
whereby two-particle processes, such as spin, charge, or 
pair excitations are taken into account by introducing 
appropriate hard-core bosons. 
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FIGURES 




FIG. 1. Visualization of the cluster perturbation approach 
for stripes: The infinite lattice is divided into unit cells. The 
unit cells differ for different stripe configurations. As shown 
here, the unit cell for a 3+1 configuration consists of two 3x6 
ladders on top of each other next to a 12-site chain. The 3-leg 
ladders on the left half of the unit-cell were diagonalized with 
a staggered magnetic field which was oriented in the opposite 
direction of the one on the right half. The hopping terms 
connecting the exactly-solved clusters (3x6 and 1 x 12) as 
well as the hoppings connecting the unit cells are included via 
the cluster perturbation technique. 



a 






b 

: / 


r 





(0,0) (JI.O) (71,71) (71/2,71/2) (0,0) 



FIG. 2. Single particle spectral function for a 2+2 
bond-centered stripe configuration at 12% doping: Compar- 
ison of CPT calculations based on exact diagonalizations of 
two 2-leg t-J ladders with J = 0.4t (a) and on exact diag- 
onalizations of two 2-leg Hubbard ladders with U = 8t (b). 
The gray scale represents the weight of spectral function at 
the specific {k, aj)-point with dark areas corresponding to high 
spectral weight. 




FIG. 3. Integrated spectral weight of bond-centered stripe 
configurations at 12% doping; comparison of CPT calcula- 
tions based on exact diagonalizations of t-J ladders (a,b; 
J = QAt) and Hubbard ladders (c,d; U — 8t): (a,c) total 
integrated weight in photoemission (n(k)), (b,d) low energy 
excitations (integrated weight in Ef — 0.2t < w < Ef)- The 
data are plotted for the whole Brillouin zone with the F-point 
in the center. The result of the stripe calculations have been 
symmetrized to account for the differently oriented stripe do- 
mains in real materials. Regions of high spectral weight cor- 
respond to white areas. 



0% 50% 0% 50% 12% 25% 12% 25% 




FIG. 4. Unit cells for different stripe configurations: (a) 
4+1 site-centered with 10 % doping; (b) 2+2 bond-centered 
with 19% doping (2x8 with 2 holes + 2x8 with 4 holes); 
(c) 2+2 bond-centered with 19% doping (2x8 with holes + 
2x8 with 6 holes); (d) 3+1 site-centered with 21 % doping. 
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FIG. 5. Electron occupation in real space before (thick 
lines) and after (thin lines) application of CPT: (a,b,c,d) ac- 
cording to Fig.(0). 




FIG. 6. Total integrated spectral weight in photoemission 
n{k): (a,b,c,d) according to Fig.([lj). 




FIG. 7. Integrated spectral weight around the Fermi sur- 
face: (a,b,c,d) according to Fig.(0). 
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FIG. 8. Density plot of the single particle spectral func- 
tion: (a,b,c,d) according to Fig.(0). 
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